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ABSTRACT 

The deflection, potential, shear and magnification of a gravitational lens following an elliptical power law mass model are investigated. 
This mass model is derived from the circular power law profile through a rescaling of the axes, similar to the case of a singular 
isothermal ellipsoid. The resulting deflection can be calculated explicitly and given in terms of the Gaussian hypergeometric function. 
Analytic expressions for the remaining lensing properties are found as well. Because the power law profile lens contains a number of 
well-known lens models as special cases, the equivalence of the new expressions with known results is checked. Finally, it is shown 
how these results naturally lead to a fast and accurate numerical scheme for computing the deflection and other lens quantities, making 
this method a useful tool for realistically modelling observed lenses. 
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1. Introduction 


Investigations in strong gravitational lensing often require an ap¬ 
proximation of the lens system by an idealised model with a 
fixed radial mass profile. In the usual formulation of gravita¬ 
tional lensing, the three-dimensional mass distribution is inte¬ 
grated along the line of sight, resulting in a two-dimensional 
surface mass density k that is responsible for the deflection of 
passing light rays (an introduction and reference can be found 
e.g. in [Schneider, Kochanek, & Wambsganssj2006j ). The initial 
model for a lens system is often a spherically symmetric mass 
distribution, which leads to a circularly symmetric surface mass 
density. Many of the classical lenses fall into this category, such 
as isothermal spheres or the Navarro, Frenk and White profile 
( jNavarro, Frenk, & White|1996l Bartelman n[l996) l. Because of 
the high level of symmetry, their lensing properties can often be 
worked out analytically. 

The next step in approximating real lens systems is to turn 
the known lenses with spherical symmetry into ellipsoids by 
rescaling an arbitrary axis of either the surface mass density 
or the deflection potential, fixing the other through the Poisson 
equation (e.g. |Kassiola & Kovner| 1993] and references therein). 
Choice of an elliptical potential - often called a pseudo-elliptical 
model - simplifies the problem, because the shear, deflection, 
and convergence can all be expressed directly in terms of deriva¬ 
tives of the potential, thus eliminating the need to solve com¬ 
plicated integrals. However, this approach can lead to unreal¬ 
istic and unphysical surface mass densities with peanut-shaped 
isodensity contours or negative values. These problems become 
more acute as the ellipticity increases, and the approach is gen¬ 
erally unsuited for axis ratios of q as 0.5 and below. 

More realistic lens models might therefore be created from 
elliptical surface mass densities. The properties of such lenses 
were described early b y|Bourassa, Kantowski, & Norton| ( |1973) > 
and Bourassa & Kantowski ( jl975[ )~ who introduced a complex 
formalism of gravitational lensing to simplify the necessary cal¬ 
culations. Fater|Schramm|(jT990!) found the equivalent deflection 
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angle in terms of two-dimensional real coordinates. However, 
due to the loss of symmetry, it is often no longer possible to find 
the lens properties of elliptical mass distributions analytically. 
Notable exceptions are the singular and nonsingular isothermal 
ellipsoids described by Kassiola & Kovner (1993)> and Kormann, 


Schneider, & Bartelmann ( |1994| > using the complex formalism, 
which feature beautifully simple closed-form solutions while 
still being widely applicable. Despite the success of this model, 
real observations often require more flexibility in the assumed 
lens profile, and no equally well-established elliptical lens with 
a simple mathematical form and numerical implementation is 
available today. 

In the present article, a general class of lenses following an 
elliptical power law profile is investigated. Previously, a number 
of numerical recipes for the nonsingular variant of this model 
have been presented ( Schneider & Weiss|l991[ Barkana|l998) l. 
In section [2] the elliptical power law profile lens is reexamined 
using the complex formalism of lensing, with the goal of find¬ 
ing compact expressions for the deflection, potential, shear, and 
magnification. In section[3] it is shown that the new expressions 
reduce to known results for a number of special cases of the 
power law profile. Finally, in section [4] it is demonstrated how 
the results naturally lead to a simple and efficient scheme for 
numerical computation. 


2. The elliptical power law profile lens 

The dimensionless surface mass density k of a lens following 
a power law profile with circular symmetry is (e.g. [Evans &| 
Wilkinson 1998 Schneider et al.|2006 and references therein) 


«(r) = 


2-t 



( 1 ) 


where 0 < t < 2 is the slope of the profile, b > 0 is the scale 
length and r > 0 is the distance from the centre of the mass 
distribution. Such a profile arises from a spherically symmetric 
three-dimensional mass distribution p(r) oc r ' 1 . The power law 
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profile lens is a versatile model; it contains as special cases the 
singular isothermal sphere for t = 1, the point mass for t = 2, 
and approximates the inner slope of the |Navarro, F renk, & White 
( |1996| > and Moore et al. |( |1999[ ) profiles for t = 0 and t = 1/2, 
respectively. 

In order to turn the profile 0 into an elliptical surface mass 
distribution, the x-axis is now stretched by a factor of q\ where 
0 < q < 1 is a constant parameter of the new profile. It is clear 
that after this operation, the formerly circular isodensity contours 
r = const have indeed become ellipses with semi-major axis r/q, 
semi-minor axis r, and axis ratio q. The resulting elliptical power 
law profile can be written as 


k{R) = 


2-t 

~Y~ 



( 2 ) 


where R is the elliptical radius defined by 

R = yjq 2 x 2 + y 2 , (3) 

i.e. the semi-minor axis of the ellipse passing through x, ij. One 
also define^ the corresponding elliptical angle 

ip = arc tan (gx, y ) (4) 

as the polar angle of x, y when transforming back to circular 
symmetry. The inverse coordinate transformation is 

x - R/q cos(ip) , 

y = R sinfv?) , (5) 


where R and ip are limited to R > 0 and ip 6 [0, 2n\, respectively. 

For the following calculation, the complex formulation of 
lensing for spheroidal mass distributions is used, as laid out by 
|Bourassa et al.| ( [T973 l and Bourassa & Kantowski|( 1975] ), with 
contributions from Bray| ( l984| l, and Schramm 1 1990!, who also 
investigated the potential of such lenses. These results were later 
used to great effect by Kassiola & Kovner ([1993 ) and even more 


Kormann et al. ( 1994|> in the analysis of the singular and 


nonsingular isothermal ellipsoid. The derivations of this section 
closely follow their outlined path. 

In complex notation, the gravitational lens equation is 


z = z~ a(z ), 


( 6 ) 


where the complex coordinate z = x + iy and complex deflection 
angle a = a x + i a tJ replace the corresponding two-dimensional 
vector quantities. Following [Schramm] ( jl990j >, one can further 
introduce the complex deflection potential i f/(z), the real part of 
which is the usual deflection potential of the real formulation. It 
is related to the complex deflection angle a(z) by the Wirtinger 
derivative (see e.g. Kormann et al.|1994 for a definition) 


v , o ^ 

a (z) = 2 — 

dz 


(7) 


where the asterisk a* denotes complex conjugation. Similarly, 
such relations can also be found for the second derivatives of the 
potential, namely the complex shear y - yi + i yi. given by 

V t da * rat 

7 (z) = -x- , (8) 

dz 

1 This definition uses the two-argument inverse tangent arc tan (x, y), 
which respects the quadrant of x, y to return the correct angle. 
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and the convergence 

da * 

^ = aF' 


(9) 


Since the latter expression must recover the surface mass density, 
it presents a direct check for the correctness of the calculations. 

The complex deflection angle for an elliptical, homoeoidal 
mass distribution was given by Bourassa & Kantow ski| ( |1975| . 
Using the elliptical radius R defined in ([3]), the deflection angle 
of an elliptical surface mass density k(R) is 


a*(z) = 2 


Vz? 


dr 

Jo 


K(r) r 


sjq 2 Z 2 - 


(1 -q 2 )r 2 


( 10 ) 


where R(z) is the semi-minor axis of the ellipse passing through 
point z. The factor in front of the integral is due to Bray| ( |T984| ; 
it ensures the correct sign of the deflection in all quadrants of the 
complex plane. Inserting the elliptical power law profile 0 into 
integral ([TO}, the complex deflection angle can be written as 


a*(z) = 


2-t b' 
q Z 
2 — t b 1 
2 q z 




rT 1/2 

z 2 ) 




1 ~q 2 R(Z) 2 , \ . 1/2 

q 2 Z 2 V 


10 

q Z 




r. Ezt\ 

2’ cf- z* ) > 


(ID 


where the change of variable r —> £ = r 2 /R(z ) 2 was used in 
the first step, and the integral representation of the Gaussian 
hypergeometric function 2 F’i(a, b; c; z) was used in the second 
step. Because its parameters are related as c = a + b + 1/2, a 
quadratic transformation of the hypergeometric function exists, 
and the complex deflection angle can further be simplified to 


a*(R,ip) 


2b 
1 + q 



1 

2 



iza 

1 +9 


e' i2 ' £ ’) , (12) 


where coordinate transformation ([5]) was applied to the complex 
variable z = x + i y. For reference, the complex deflection angle 
before complex conjugation is 


a(R, <p) 


2b 
1 + q 



1 

2 


1 ~g -i2y>' 
1 +q 


(13) 


This result is a beautiful factorisation of a into its (elliptical) 
radial and angular parts; a fact that will be exploited in section[4] 
to calculate the deflection quickly, in spite of the hypergeometric 
function it contains. 

Taking inspiration from the singular isothermal ellipsoid, 
one finds that a potential which solves equation 0 for complex 
deflection © is given by 



za'fz) + z* a(z) 
2 


(14) 


In real coordinates, this expression reduces to the more familiar 


xa x + ya u 

<A(x,y)= . ; 


(15) 


It is clear that this potential can be calculated at very low cost 
if the deflection angle a has already been found. Checking the 
result is tedious; one possibility is to substitute e -1 ^ = (qx-iy)/R 
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Fig. 1. Isocontours of the deflection potential ^ (black) and surface mass density k (grey) for an elliptical power law profile lens in physical 
coordinates x, y. For the purpose of illustration, the contour levels are equally spaced along the diagonals. The slope of the power law profile varies 
from t = 0.25 (left) to t = 1.75 (right) in steps of 0.5. The axis ratio varies from q = 0.8 (top) to q = 0.2 (bottom) in steps of 0.2. 


and e~ l2t< ’ = (qx - i y)/(qx + i y) in deflection ( [T2 | i and express 
the Wirtinger derivative in ([7]) in terms of real coordinates as 
d/dz = (d/dx - id/dy)/2. 


The complex potential ( fl4| evidently has no imaginary com¬ 
ponent. In its form ( p~5| ), it can therefore be taken to be the usual 
real deflection potential satisfying a = Vi/r (|Schramm||l990|l. 
Fig-0 shows the potential for various settings of the power law 
slope t and axis ratio q. It has been noted that the potential is al¬ 
ways "rounder", i.e. less eccentric, than the surface mass density, 
and this effect is clearly visible. 


Calculating the shear from equation ([8]) is straightforward, 
and results in 

7 * Ct*(z) 

y\z) = -K(z)- +(l-f)— . (16) 

z z 

For computations, the shear is more clearly expressed in terms 
of the physical polar coordinates r and 8 as 

y = -e i2s k + (\ - t)t' e — , (17) 

r 

once more recalling the singular isothermal ellipsoid (f = 1). As 
in the case of the potential, the shear is readily calculated once 
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the surface mass density k and complex deflection angle a are 
known. 

Given shear y and convergence k, it is now possible to find 
the (inverse) magnification — (1 - a :) 2 - \y\ 2 of the elliptical 
power law profile lens, which is 


i / x a x + u a u \ 

A 7 = l-2/t(l-(l-0 x r / y J-(l-f)- 


\a\ 2 


(18) 


Here r denotes the radius r 2 = x 1 + y 2 in physical coordinates. 
Using expression it is possible to determine the critical 
lines p 1 = 0 and corresponding caustics of the lens. Because the 
resulting equations contain the hypergeometric function, their 
solutions are found numerically. The critical lines and caustics 
for a number of settings of the power law slope t and axis ratio q 
are shown in Fig. [2] 


3. Special cases 

In order to check the presented solutions for the deflection ( fl3} , 
potential < p~4] ), shear and magnification ( p~8| ), it is useful to 
compare them to known special cases of the power law slope t. 

The singular isothermal ellipsoid is a power law profile lens 
with slope t = 1. The hypergeometric function can be simplified 
in this case, yielding the complex deflection angle in physical 
coordinates 


a(x, y) = 


2b 

V 1 -q 2 


arctan 


Vi ~ q ^qx + i y 
y/l + q \jqx - i y t 


(19) 


Taking the real and imaginary part and using the sum formula 
for the inverse tangent reduces the expression to the commonly 
used form as reported e.g. by Schneider et al. ( 2006} . It is trivial 
to check the potential <// = xa x + ya y , shear y = -kz* /z, and 
magnification p~ l = 1 - 2k. 

The limit t —> 2 takes the power law profile to a point mass. 
The hypergeometric function exists in the limit, resulting in the 
complex deflection angle 


a{x, y) = 


b 2 (x + i y) 
qr 2 


( 20 ) 


As expected for a point mass, the deflection is always circularly 
symmetric (note r instead of R), and the axis ratio q appears only 
in the form of a constant of normalisation, making the Einstein 
radius of the lens r E = bq l/2 . While the potential ([14} cannot 
be evaluated in the limit t —» 2, both the shear y = -It/( qr 2 ) e l2e 
and the magnification p 1 = 1 -b 4 /(q 2 r 4 ) are easily checked and 
agree with the known results for the corrected Einstein radius. 

The final case of interest is the uniform critical mass sheet in 
the limit t —» 0. In this case, the deflection becomes 


a(R, <p) 


1 + q 


Re 1 ' 


( 21 ) 


which has a clearly unphysical dependency on ip, as the surface 
mass density k = 1 has circular symmetry. However, the same 
result ( [2T] > is obtained if k = 1 is inserted into deflection ( [TO} 
from |Bourassa & Kantowski ( |1975) ) directly. The problematic 
result is therefore a limitation of the formalism for elliptical mass 
distributions, and not the elliptical power law profile lens. 
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4. Numerical evaluation 


The deflection angle ( fTJ} is at the heart of the elliptical power 
law profile lens, as the potential ( fi~4| and shear ( | 1 6} can all be 
expressed in terms of it. It is therefore necessary to find a fast 
method for evaluating the contained hypergeometric function to 
make this a useful model for numerical simulation and modeling 
purposes. For this, it is useful to separate the radial and angular 
parts of the complex deflection angle (fl3}, which becomes 


a(R, ip) 


2b 
1 + q 



( 22 ) 


The angular dependency of a is contained in the function 

w(q>) - z lip iF\ (l, |;2 — —f e i2(P ) , (23) 

where / = is the second flattening of an ellipse with axis 
ratio q. Because 0 < q < 1, the range of / is 0 < / < 1, and it 
follows that the hypergeometric function in expression ([23} has 
a convergent series representation 


«(</>) = 'Yj 


n =0 


r (2-|)T(n + |) 
r (^)T(n+2 - {) 


(-/)" 


£ i(2«+l)(o 


(24) 


This, on the other hand, is nothing but a Fourier-type series 


<4<p) = Y a >^ ii2n+1 ^ 

n=0 


(25) 


containing only positive and odd terms 2n + 1, the coefficients of 
which are given by 


a„ 


F(2-I)r(n + 0 

r(i)r(n + 2-|) 


(26) 


Since this is a hypergeometric series, the ratio of two subsequent 
series coefficients is simple. 


a„ _ 2 n - (2 - t) 

a „-1 2 n + (2 — f) 


(27) 


It is clear that the magnitude of the series terms drops off almost 
geometrically, with an asymptotic rate of /. This behaviour is 
shown in Fig.[3]for various settings of the power law slope t and 
axis ratio q. 

The ratio ( [27] ) can be used to iteratively calculate the terms 
of series ([25}. Introducing symbols for the summands as 

oo 

oj= Y Q " ’ ( 28 > 

n =0 


the /j-th term Q„ is related to the previous one as 

= 2n - (2 - t) i 2l p q (29) 

J 2n + (2 - 0 

This reduces the calculation of the deflection to summation and 
(complex) multiplication. Furthermore, since the computation is 
iterative, it can be continued easily until the desired precision or 
accuracy in the deflection is reached. 

Instead of the numerical scheme given by equations ( [28} 
and ([29}, it can be advantageous to forego the use of complex 
numbers, particularly in computer implementations. Just as the 
two-dimensional deflection angle a is given by the real and 
imaginary components of the complex deflection angle a , the 
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Fig. 2. Critical lines (black) and caustics (grey) of the elliptical power law profile lens. Also shown is an ellipse with semi-minor axis equal to the 
scale length b and axis ratio q (red, dotted). For the purpose of illustration, the scale length b decreases as q >/2 . The slope of the power law profile 
varies from t = 0.25 (left) to t = 1.75 (right) in steps of 0.5. The axis ratio varies from q = 0.8 (top) to q = 0.2 (bottom) in steps of 0.2. 


complex function to for the angular dependency corresponds to 
a two-dimensional vector a> in the real formulation. Expanding 
the exponential in Fourier series ( [25] ) into its real and imaginary 
parts, the components of a> can be written as the individual series 


oo 

= 2 a n cos ((2 n + 1 )ip) , 

n =0 

oo 

^ a„ sin ((2 n + 1) ip) . (30) 

n =0 


Just as in the complex case, the real components of the deflection 
angle can be calculated iteratively. Writing 


oo 



n=0 


oo 

oj y = ^Q^„, (31) 

n-0 

the relations for the n -th terms can be found by understanding 
the complex product e l2t Q„ i in expression ( [29] ) as the two- 
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Fig. 3. Absolute value of the series coefficients a„ for the elliptical 
power law profile lens. Shown are different values t = 0.5 (dotted), 
f = 1.0 (solid), and t = 1.5 (dashed) for the slope of the power law and 
axis ratios from q - 0.8 (bottom) to q = 0.2 (top) in steps of 0.2. 


for the quick numerical evaluation of the deflection angle, and 
consequently other lensing quantities. 
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dimensional matrix multiplication 
2 n - (2 - t) 

«« = -/ - J. ’ R^O-i , (32) 

2 n + (2 — t) 

where is the rotation matrix. In components, this is 

Qx.n = -f + (7 _ f) ( C0S ( 2 V) Q -vi-i - Sin(2^)0 !/i „_ 1 ) , 

&y,n = ~f 2n + ( 2 - f ) ( Sln( ^) + COS(2 ip) , ) . (33) 

This prescription can make evaluation of the deflection suitably 
fast to use the elliptical power law profile lens for ray-tracing in 
model fitting and Monte Carlo inference, where a large number 
of calculations are necessary. 


5. Conclusion 

It was shown how the power law profile (|T]> for lensing can be 
turned into a surface mass density Q with elliptical isodensity 


contours. Following the groundbreaking work by Bourassa et al. 
(|1973|> and Bourassa & Kantowski] (|T975]> on ellipsoidal lenses, 


and Kormann et al. ] 1994] > on the singular isothermal ellipsoid, 
the properties of the elliptical power law profile lens have been 
derived using the complex formalism of gravitational lensing. It 
was found that the deflection angle ( fT3] l can only be expressed 
in terms of the Gaussian hypergeometric function nf) (a, b; c; z), 
thereby severely limiting the analytical tractability of the lens 
equation. Further results for the potential ( p~4] >, the shear ( f]~6l >, 
and the magnification ( fl8| > could all be given in terms of the 
complex deflection angle, which therefore becomes the central 
quantity for computations. These results have been checked in 
section[3]using a number of known special and limiting cases of 
the power law profile. Finally, a recipe was given in section [4] 
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